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ABSTRACT 



The discovery of the 7-ray pulsar PSR J1836+5925, powering the formerly uniden- 
tified EGRET source 3EG J1835+5918, was one of the early accomplishments of the 
Fermi Large Area Telescope (LAT). Sitting 25° off the Galactic plane, PSR J1836+5925 
is a 173 ms pulsar with a characteristic age of 1.8 million years, a spindown luminosity 
of l.lxl0 34 erg s _1 , and a large off-peak emission component, making it quite unusual 
among the known 7-ray pulsar population. We present an analysis of one year of LAT 
data, including an updated timing solution, detailed spectral results and a long-term 
light curve showing no indication of variability. No evidence for a surrounding pulsar 
wind nebula is seen and the spectral characteristics of the off-peak emission indicate 
it is likely magnetospheric. Analysis of recent XMM-Newton observations of the X-ray 
counterpart yields a detailed characterization of its spectrum, which, like Geminga, is 
consistent with that of a neutron star showing evidence for both magnetospheric and 
thermal emission. 

Subject headings: gamma rays: general; pulsars: general; pulsars: individual (PSR 
J1836+5925) 



Introduction 



Since its discovery by EGRET (ILin et al 1 ll992h . the bright high-energy 7-ray source GRO J1837+5J 
defied straightforward identification. It was reported as a persistent source with a varying flux of 
3-8 x l(T 7 ph cm s and a relatively hard spectrum of photon index 1.7 in non-consecutive, 
typically 2-3 week-long, observing periods. Its location at hi gh Galactic latitude in coni unc- 



tion with early reports of 7-ray variability (later questioned by iNolan et al.l 1 19961 : iReimer et al 



200 ll ) suggested it might be a blazar. However, the lack of a radio-bright counterpart, common to 
EGRET-detected blazars, cast doubts on such an interpretation. 



With the detection of faint X-ray counterpart candidates in the error contour of 3EG J1835+5918 (jReimer et al 



2000), the interpretation focused increasingly on a nearby radio-quiet neutron star. The complete 



charac teriz ation of all but one of the ROSAT HRI X-ray sources was presented by iMirabal et al 



(|2000j) and lReimer et al.1 (|200ll ). who singled out RX J1836. 2+5925 as the most probable counter- 
part of 3EG J1835+ 5918. Subaru/FQC AS observations in the B- and U-bands proposed possible 
optical counterparts (ITotani et al.1 120021). w hile Hubble Space Telescope observations set an optical 
upper limit of V > 28.5 (jHalpern et all2002 ) . A scenario of a thermally emitting neutron star which 
was either older or more distant than the archetypal radio-quiet 7-ray pulsar Geminga emerged 
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as the most plausible explanati on for the source ( Halpern &: Rudermanlll993l ; iBignami fc Caraveo 



19961 ; iMirabal &: Halpern! 120011). with an upper limit on the distance of 800 pc, determined from 
X-ray observa t ions (IHalpern et al.l 120021 ) . Using Chandra observations separated by three years, 



Halpern et al.l (|2007l ) were also able to determine an upper limit on the proper motion of 0.14"per 
year, or v t < 530 km s" 1 at 800 pc. However, a ti ming signature, which would settle the n ature 



of this source, was never found in the EGRET data (IChandler et al. 

boogj. 



in repeated observations by Chandra (Halpern et al 



NRAO's Green Bank Teles cope (GBT) (IHalpern et al.ll2007f ). Upper limits from Very High Energy 



2001 



Ziegler et al.ll2008h . nor 



2007), nor in a 24- hr observation with 



(VHE) 7-ray observations ( Fegan Sz Weekes 20051 ) determined that the peak of emission must be 



at GeV gamma rays. Recently, AGILE reported marginal flux variability in their 2007-2008 data, 
arising from several non-detections in a period of long uninterrupted c overage, along with a flux 
level significantly lower than what was previously reported by EGRET (jBulgarelli et al.l 12008). 



3EG J1835+5918 was a target of pointed observations during the 60-day commissioning pe- 
riod prior to the start of normal science operations of the Large Area Telescope (LAT) aboard 
the Fermi Gamma-ray Space Telescope. During these observations, the LAT accumulated pho- 
tons from this source at a rate approximately twice as high as during regular survey-mode oper- 
ations, thus facilitating the detection of 7-ray pulsations. The discovery and initial ti ming of the 



pulsar, PSR J1836+5925, using the first five months of LAT data, were reported in lAbdo et al 



(|2009al ). Here we present the phenomenology emerging from one year of LAT observations of 
PSR J1836+5925, including energy-dependent pulse profiles and phase-resolved spectroscopy. 



2. Gamma-ray observations and data analysis 



The LAT is a pair conversion telescope, sensitive to gamma rays with energies from 20MeV 
to >300 GeV. Gamma rays in the LAT are recorded with an accuracy of < 1 fis. The LAT has an 
on-axis effective area of 8000 cm 2 , a field of view of ~2.4sr, and an angular resolution of ~0.8° 68% 
containment at 1 GeV (jAtwood et al.ll2009i ). 



2.1. Timing analysis 

We have derived a precise timing solution of PSR J1836+5925 using data from 2008 June 30 
to 2009 June 30 (MJD 54647.4-55013.0). We selected photons with E > 170 MeV offset from the 
source direction by no more than 1.6°, a radius chosen to maximize the pulsed significance, and 
used the LAT Science T00H gtbary in its geocenter mode to correct the arrival times to terrestrial 
time (TT) at the geocenter. We generated a total of 22 pulse times of arrival (TOAs), each covering 
roughly two weeks of data, and obtained pulse profiles by folding the photon times according to 



available at 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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a provisional ephemeris using polynomial coefficients generated by Tempo2 (jHobbs et al.l 12006) 
in its predictive mode (assuming a fictitious observatory at the geocenter). The TOAs were then 
measured by cross correlating each pulse profil e with a template consisting of two gaussians, de- 



rived from the data set above (|Rav et al 



2010j). The timing model, fit using Tempo2, included 



position, frequency (v), and frequency derivative (z>). Table [T] lists the results of our fit. The 1.3 ms 
RMS residual to the fit is comparable to the mean TOA measurement uncertainty of 1.2 ms and 
significantly smaller than the 5.4 ms resolution of our 32 bin light curve. The reduced x 2 °f our 
timing fit is 1.9. With v = 5.77 Hz and v = —5 x 10~ 14 Hzs -1 , we derive a characteristic age of 1.8 
million years and a spindown luminosity of 10 34 ergs -1 . Our best fit location is RA=18:36:13.75(3), 
Dec=+59:25:30.3(6), which is 0.35" from RX J1836. 2+5925, well within the statistical uncertainty 
of the timing fit, securing the association between PSR J1836+5925 and RX J1836. 2+5925. 



2.2. Light curve 

We explored the pulsar light curve in different energy bands by selecting events with energies 
above 100 MeV from an energy-dependent region of interest (ROI), defined as 6 = 3.4 o (£ , /100MeV) _ ' 
with a minimum (maximum) radius of 0.35° (2.1°). The rotation phase of each event is calculated 
using the truncated Taylor series expansion: <p = </>o + v{t — Tq) + \v{t — To) 2 , where To is the 
reference epoch of MJD 54800 and (f>o is the reference phase at To, which we define as </>o=0.55. 
Figure Q] (top panel) shows the folded light curve of the pulsar for energies above 100 MeV. The 
light curve has two distinct peaks and is well fit by a constant plus two gaussians centered at phases 
(j) = 0.26 and (j) = 0.77, with their means separated by 0.51 ± 0.01 in phase. The "pulsed fraction" 
(determined by integrating the contribution from the two gaussians) is energy dependent, as can 
be appreciated from Figure [TJ It has an unusually low value of ~ 26 ± 2 % for energies above 100 
MeV, with the remaining >70% coming from the constant term. While the pulsed fraction does 
increase at higher energies, up to ~45% above 1.5 GeV, it is always less than 50%, regardless of the 
cuts chosen. 

We identify the following intervals: first peak (FP): 0.105< <p <0.405, second peak (SP): 
0.632< <p <0.904, bridge (BR): 0.459< <j) <0.597, and off-peak (OP): 0.938< <j> <0.053. The lower 
panels show the folded light curve in 5 different energy intervals: 0.1-0.3 GeV, 0.3-1 GeV, 1-3 GeV, 
>3GeV, and > 5 GeV (dark histogram in the second panel). 



2.3. Spectral Analysis 

We performed the spectral analysis using data collected during the sky survey: 2008 August 4 to 
2009 June 30 (MJD 54682-55013). While the LAT data taken during the commissioning period are 
adequate for timing analyses, several of the configuration settings may have had a modest effect on 
the energy resolution and reconstruction, so we exclude these data for the spectral analysis. We also 
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exclude events with zenith angles greater than 105° to minimize the contamination from gamma 
rays from the Earth's atmosphere. A phase-averaged spectrum was obtained with an unbinned 
maximum likelihood analysis, using the LAT Science Tool gtlike (jAbdo et al.l l2009bl ) and the 
"Pass 6 v3" instrument response functions (IRFs). We used energies > 200 MeV. The diffuse 
emission from the Milky Way was modeled using gll_iem_v02 1 while the isotropic extragalactic 
diffuse emission and residual instrumental particle backgrounds were modeled together using the 
currently recommended "template spectrum" isotropic_iem_v02 1 . We extracted photons from a 
15° radius ROI centered on the coordinates of RX J1836. 2+5925, in order to properly account for 
the contributio ns of other gamm a-ray sources in the vicinity. All sources from the LAT first-year 
source catalog (jAbdo et al.ll2010l ) included in our ROI were fit with a simple power law, while for 
the pulsar itself we used a power law with an exponential cutoff: 



diV 
dE 



TT7 = ^GeV ex P 



E 



E, 



cutoff 



(1) 



where T is the photon index, E cuto g the cutoff energy, and K the normalization (in units of ph cm~ 2 
s _1 MeV -1 ). A super-exponential cutoff (n = 2) was ruled out, compared to a simple exponential 
(n = 1) cutoff, at 8a significance so we set n equal to 1 in further analysis. Note that physical 
motivation for the use of a power-law spectrum with an exponential cutoff is underpinned by this 
being approximately the form expected for curvature or syn chrotron radiation from both monoen- 
ergetic electrons (e.g. see Eq. (24) of lHarding et al.ll2008l ). and electrons with a distribution of 
Lorentz factors up to some maximum value. A superexponential cutoff, on the other hand, would 
be expected in polar cap m odels, due to single photon pair production attenuation in strong mag- 
netic fields near the surface (jNel k, de Jagerlll995l : lDaugherty fe Hardinglll996l : lRazzano Harding 



20071 ). The results of our spectral fits are summarized in Table [2j The quoted errors are statistical 



only. The effect of the systematic uncertainties in the effective area on the spectral parameters 
is 5T = (+0.3, -01), <5£ r , lt off = (+20%, -10%), 5Fi 00 = (+30%, -10%), and 5Gi 00 = (+20%, 
-10%) (|Abdo et al.l 12009^ . 

Next, we studied the energy spectrum in the various phase intervals defined in Section \2.2[ The 
spectral parameters of the pulsar were allowed to be free, while those of the other sources were fixed 
to the values obtained in the phase-averaged analysis. Figure [2] shows the four spectra and Table [2] 
summarizes the results, normalized for the different phase intervals, pointing to mild variations of 
the photon index T and cutoff energy 22 C utoff over the four phase intervals, with the off-peak region 
characterized by a softer spectrum. We also carried out a phase-resolved spectral analysis by taking 
15 equal-counts bins of ~650 photons. As in the previous analysis, the parameters of all the sources 
in the ROI were fixed at the values obtained in the phase- averaged analysis. Figure [3] shows the 
results of our analysis. The top panel illustrates the evolution of the cutoff energy, while the bottom 
panel shows the change in photon index with phase. Insufficient photon statistics prevent us from 
investigating the apparent spectral changes within the peaks in any finer detail. 
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2.4. Variability analysis 

2.4-1- Pulse profile variability 

We checked for variability in the pulse profile shape using non parametric tests. First, the 
data were divided into time segments with equal counts (from 2 up to 32 segments). For every pair 
of segments, light curves were compared against the null hypothesis that both were drawn from 
the same parent distribution by performing Kolmogorov-Smirnov tests. There was no instance 
in which the pairwise comparison resulted in the null hypothesis being rejected even at the 80% 
confidence level, x 2 tests yielded the same results, resulting in an overall gaussian distribution for 
the normalized residuals. We repeated the tests for light curves with 8, 10, 12, 15, 20, and 30 bins, 
confirming in every case that there is no hint of variability on timescales longer than a week. 



2-4-2. Flux variability 



In order to check the long term stability of the source, we computed the 0.1-100 GeV photon 
flux in 5-day time bins using the LAT Science Tool gtlike. First, we constructed a sp ectral model 
file including all the LAT first- year catalog point sources catalog (jAbdo et al.l |2010| ) in our ROI, 
the Galactic diffuse emission and the isotropic background. Since AGN in the region of interest 
can be variable, the spectral parameters of all point sources were set free while the Galactic and 
isotropic background were fixed to the values obtained in the first phase- averaged analysis. After 
running gtlike, we checked the fit result for sources with large uncertainties. After removing these 
sources with low significance (< 2a), we ran gtlike again to obtain the flux of PSR J1836+5925 for 
that 5-day bin. Figure 0] shows the re sulting fluxes an d stati stical uncertainties. We used five-day 
bins as this was the interval chosen by iBulgarelli et al.l (|2008l ). facilitating the comparison with the 
AGILE results. Assumi ng a constant flux, a y 2 tes t gives a value of 66.2 with 65 degrees of freedom. 
The variability index of iMcLaughlin et al.l (|1996l ) is V = 0.37, consistent with no variation. The 
weighted standard deviation of the flux is 16% of the average. If the flux is assumed to change 
linearly with time, the slope is consistent with zero, with a 68% limit of 6.8% change from beginning 
to end of the data. Fitting a sinusoidal variation produces an amplitude of (2.5 ±3.9)% of the mean 
flux, consistent with zero. We note t hat while the overall flux seen by the LAT agrees with that 



reported by the EGRET experiment (IHartman et al.lll999), it i s somewhat higher than what has 



been reported by the AGILE experiment (IBulgarelli et al 



2008 



Pittori et al J 12003 ). 
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3. Multi-wavelength Analysis 
3.1. Radio search 

We used the LAT ephemeris to search anew for radio pulsations from PSR J183 6+5925. We 



folde d the 24 hr dataset obtained at the GBT in 2002 December (for details, see lHalpern et al 



20071 ) modulo th e predicted per iod, while searching in dispersion mea sure up to DM = 100 cm 3 pc, 



using PRESTO (|Ransomll200ll l. For the distance range 250-800 pc (lHalpern et al]l2007T l. the DM 
predicted by the NE2001 electron density model (jCordes &: Lazioll2002l ) is 2-9cm _3 pc. Although 
there is no evidence in the LAT 2008-2009 timing solution for rotational instabilities, and we do 
not expect a large degree of timing noise from such a relatively old pulsar, we also did a small 
search in period about the nominal value. No radio pulsations were detected. For an assumed 
pulsar duty cycle of 10%, our long observation at a frequency of 0.8 GH z yields a flux densit y of 
Sq.8 < 7/xJy (this is a significant improvement over the limit presented in lHalpern et al.l (|2007l ) for 
the same data because we are now searching for a known period, allowing for a lower signal-to- 
noise ratio detection threshold). Converted to the more usual pulsar search frequency of 1.4 GHz 
with a typical spectral index of -1.6, S1.4 < 3 /iJy. The implied luminosity is L1.4 = Si^d 2 < 
0.002dQ 8 mJy kpc 2 . This is at least an order of magnitude smaller than the least luminous radi o 
pulsar, PSR J1741— 2054, originally discovered in gamma rays by the LAT (see lCamilo et al.ll2009l ). 
suggesting that if PSR J1836+5925 is an active radio pulsar its beam probably does not intersect the 
Earth. One caveat to this conclusion is that scintillation caused by the interstellar medium could be 
quite significant for this observation of this pulsar: depending on its actual DM within the expected 
range, the characteristic scintillation bandwidth and timescale at 0.8 GHz predicted by the NE2001 
model may be greater than the observation bandwidth and time. If so, the received flux density 
of the pulsar during the observation may not reflect its intrinsic average. To address this potential 
concern, we did one extra observation with the GBT. On 2009 October 24 we recorded data from 
a bandwidth of 100 MHz centered on 350 MHz for 2.0 hr using GUPpfl. A gain, no pulsations were 
detected from PSR J1836+5925. For the same assumed duty cycle, the flux limit was 55 /iJy With 
the same assumed spectral index, this corresponds to 14 fiJy at 0.8 GHz. While this recent 350 MHz 
observation was thus nominally only half as sensitive as the earlier 24 hr observation at 820 MHz, it 
was still an extremely deep observation, and much more immune to scintillation effects, rendering 
our earlier conclusion valid: for all practical purposes, PSR J1836+5925 is a "radio quiet" pulsar. 



3.2. X-ray observations 

We studied the X-ray counterpart of PSR J1836+5925 using two X MM-Newton observa tions 

and 



taken on 2008 M ay 18 and 2008 Jun e 25 (15 ks each). Both the EPIC/pn (IStriider et al. 



2001 



the EPIC/MOS (ITurner et al.ll200ll ) cameras were operated in their Full Frame mode, using the 



2 https://wikio.nrao.edu/bin/view/CICADA/GUPPiUsersGuide 



- 10 - 



thin optical filter. We focused on spectroscopy (the time resolution is not adequate to search for 
pulsations). Data reduction and analysis were performed with the XMM-Newton Science Analysis 
Software (SASv8.0). Owing to the lack of variability between the two epochs, the two data sets 
were merged, resulting in 24.4 ks, 30.7 ks and 30.8 ks of good exposure in the pn, MOS1 and MOS2 
cameras, respectively. The background-subtracted 0.2-3 keV count rate of the source, as extracted 
from a 15" radius circle, is 0.023 ± 0.001 cts s -1 , 0.0043 ± 0.0004 cts s _1 , and 0.0040 ± 0.0004 cts 
s _1 in the pn, MOS1 and MOS2 cameras, respectively. Background accounts for ~ 30% additional 
counts. 

We performed simultaneous fits to the pn, MOS1 and MOS2 spectra using the XSPEC vl2.4 
software. The X-ray spectrum cannot be described by a pure blackbody model (xt = 3.70, 40 
d.o.f.). A simple power law is possibly consistent with the data (xt = 1-39, 40 d.o.f.), however, the 
best fit requires a rather large photon index (r = 3.0 ± 0.2) and a very low Njj of < 2 x 10 19 cm -2 
(errors are at 90% confidence level for a single parameter). 

The combination of a blackbody and a power law yields a better fit (xt = 0.71, 38 d.o.f.). 
The best fit model features an absorbing column Njj < 2.7 x 10 20 cm -2 (the best fit value is 0), a 
blackbody temperatur^fl kT = 59l^ 7 eV, and a power law photon index T = 1.7 ± 0.3. The total 
observed flux in 0.2-5 keV is 5.5 x 10 -14 erg cm~ 2 s _1 . Table [3] summarizes the results of our best 
spectral fit. 



While our flux values are not dissimilar from those of lHalpern et al.l (|2002l ). the black-body 
temperature, as well as the Njj values can now be constrained directly on the basis of the XMM 
data. 

The rather high value of the blackbody temperature, coupled with the very low Nh (consistent 
with 0) point to a small emitting surface at a relatively low distance. For the best fitting tempera- 
ture, an emitting surface of 1 km radius would imply a 450 pc distance which would scale to 300 pc 
for a 50 eV temperature. In no way can emission from the entire neutron star be invoked since it 
would imply a distance in excess of 3 kpc, not compatible with our very low Njj. In such a scenario 
the contribution of the thermal emission from the bulk of the surface should be negligible within 
the EPIC band. This requires^ a surface temperature lower than ~ 25 eV (~ 30 eV), assuming a 
10 km emitting radius at 300 pc (450 pc), which is consistent with expectations for a ~ 10 6 yr old 
neutron star. Such a cooler thermal component would also be consistent with the deep HST upper 
limits. 

The overall similarity of the X-ray spectrum of PSR J1836+5925 and Geminga is apparent 
in Figure [5j where both the XMM spectra and curves fitting the phase-averaged LAT spectra 



3 We quote blackbody temperatures and emitting radii as measured by a distant observer throughout the paper. 

4 In order to get a rough estimate, we fixed all spectral parameters of the blackbody plus powerlaw model to their 
best fitting values, leaving Nh as a free parameter (with a maximum allowed value of 2.7 x 10 20 cm~ 2 ), and we added 
a second blackbody component to account for surface emission. 
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for both these sources are depicted. This similarity applies to both the X-ray and 7-ray spectra 
individually, to the LAT-band turnovers, to the offset between their fluxes in each waveband, and 
therefore to the oy erall multiwavelength impression (for the detailed LAT results on Geminga, see 



Abdo et al.l 12009a ). This broadband picture clearly illustrates that PSR J1836+5925 resembles 



Geminga in its high energy components, a character that will guide future spectral modeling. The 
extrapolation of the XMM power law tails up to the LAT band highlights the disparity between 
the X-ray non-thermal indices T ~ 1.7 and the LAT band power law indices T ~ 1.3. This property 
suggests that some as yet undetectable spectral structure or feature must exist in the 20 keV - 
100 MeV band, perhaps due to a transition between components spawned by different radiative 
processes. The structure may be a simple flattening, or something more complex, however it 
should be unlike the steepeni ng seen in the broadband X-ray/7-ray spectrum of the younger Vela 
pulsar (jStrickman et al.lll999l ). The diagnostic potential enabled by the detection of such spectral 
structure motivates the development of future sensitive spectroscopy telescopes in the hard X-ray 
and soft 7-ray bands. 

We also re-analyzed archival Chandra data taken in high time resolution, to sear ch for possible 



X-ray pulsations. Using 118 ks of HRC-S data (the same dataset described by lHalpern et al 
20071 ). we extracted ~ 790 source counts. No significant modulation is apparent in a 10 bin phase 



histogram folding the events with the extrapolated LAT timing solution. Although we expect no 
significant timing noise from this pulsar, a search for pulsations was also performed around a narrow 
range of the expected period (0.1732-0.1733 s), but no significant signal was detected. Following 
Vaughan et al.1 (119941 ) . we set an upper limit of 40% on the pulsed fraction (at 99% confidence 
level), assuming a sinusoidal modulation. 



4. Discussion 

The discovery by the LAT of PSR J1836+5925 confirmed the long-held suspicion that 3EG J1835+5918 
was a nearby Geminga- like pulsar. Its characteristic age of 1.8 million years agrees with expecta- 
tions that it should be significantly older than Geminga, given the soft X-ray spectrum of its X-ray 
counter part, the absence of o ptical and radio emission, and the measured upper limits on its proper 



motion (jHalpern et al.l 120071 ). We see no evidence for time variability of either the source flux or 
the pulse profile shape over the 11 months of observations. Our measured pulsations explain why 
this pulsar proved rather difficult to detect: the small (~30%) pulse fraction and relatively large 
duty cycle made blind searches of EGRET data futile. Indeed, it necessitated the LAT pointed 
observations to get an early pulse detection. Our detection, in turn, raises several puzzles - why 
does this object have such a large off-peak component and how does the relatively low spindown 
power produce an apparently large 7-ray luminosity for even modest distances? 

While other p ulsars show detectable emission throughout most of the pulsar period (e.g. see 
Abdo et al.ll2009bl ldr). that of PSR J 1836+5925 is particularly bright, and provides a good oppor- 



tunity to test the nature of these off-peak components. Since the source at pulse minimum is 
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unresolved, and since the pulse minimum spectrum demands a 2-3 GeV cutoff, we conclude that 
we are not seeing evidence for a surrounding pulsar wind nebula (PWN). The lack of extended 
emission in the deep Chandra images also implies that there is no bright PWN. Thus the off-peak 
flux is likely magnetospheric. Since the T ~ 1.6 off-peak spectral index is substantially softer than 
that of the rest of the profile, we tested whether a second, spatially unresolved pulsar could con- 
tribute the off-peak flux. We searched for puls ations from the source by applying the standard 
time-differencing technique (lAtwood et al.ll2006l ). masking the frequency of PSR J1836+5925. We 
used a maximum frequency of 64 Hz, and a long time-difference window of ~12 days. We found 
no evidence for pulsations at any other frequency. The precise sensitivity of the blind search is 
still not completely un derstood, however, a comparison of the blind search pulsars discovered so 
far (jAbdo et al.ll2009al ) and the known radio pulsars detected by the LAT, suggests that the blind 
search is approxi mately 2-3 times less sensitive than a standard pulsation search using the known 
timing solution ( Abdo et al.ll2009fl ). This results in a 5a limit on the pulsed flux of ~ 2 x 10~ 7 
cm -2 s _1 for another putative pulsar at this location. We conclude that PSR J1836+5925 emits 
over half its flux in a nearly constant component with an exponentially cut-off spectrum which is 
softer than the peaks of the profile. 



The distance inferred from the observed 7-ray flux, F 7 = 6 x 10 



-in —2 - 

erg cm s 



depends on 



the intrinsic luminosity (L 7 ), the beam geometry, and the line of sight along which we sample the 
anisotropic emission. To account for anisotropy, we parameterize the relation between the observed 
flux and the true luminosity by the "flux conversion factor" = L~ / /4ird 2 F~ / , whose estimation we 
discuss in the next paragraph. The intrinsic luminosity can be inferred from the spindown luminos- 
ity of the pulsar, if we know the efficiency 77 = L^/E. We therefore have d = (E j 'A-kF^) 1 / 2 (77 / fn) 1 ^ 2 . 
It has been argued that the efficiency of 7 -ray emission, and the fraction of the op en zone participat- 
ing in the gaps, gro ws with decreasing E (jRuderman fc Chengll988l : lAronslll996l ). and observations 
support 77 cx E~ 1 / 2 (jAbdo et al.ll2009fl ). We adopt r\ = C (£^/10 33 ergs _1 ) _1 / 2 , wher e C is a slowly 



varyi ng function of order unity which depends on the details of the physical model (jWatters et al 



2009). The observed E = 1.1 x 10 34 erg s 1 then implies an efficiency of 77 ~ 0.30. Using this 



efficiency, along with the known values for 7-ray flux F~ and spindown luminosity E, our estimate 

1 /o 

for the pulsar distance becomes d 



215/ Q 1/2 pc. 



The factor /q depends sensitively on the emission model, on the inclination of the pulsar spin 
axis to the line of sight (£) , a nd on the inclination of the magnetic pole with respect to the spin axis 
(q). Models are described in lWatters et al.l (|2009l ); for the "Two Pole Caustic" (TPC) model, pulse 
separations A = 0.5 occur in two r egions: a > 85°, C 5 -, 60° or a < 60°, £ > 85° (near the axes 
in the magenta zone of Figure 3 in IWatters et al.l l2009i ) . The former solutions are, however, not 
satisfactory as they have weak bridge fluxes, at least for models with thi n radiating surfaces. Thicker 
emission zones can produce additional bridge flux ([Venter et al.ll2009i ). The large Q solutions can 
indeed have substantial off-peak flux arising at modest r < 0.2 r^c altitudes (where r^c = cP/2it 
is the speed of light cylinder), especially for efficiencies 77 < 0.2. For the outer gap (OG) model, 
only a few £ > 80°, a < 30° models give the observed A for highly efficient pulsars (77 ~ 0.2). 
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These have relatively large off-peak fluxes, arising from large r > 0.5 r^c altitudes. 

While both models have acceptable large Q solutions, for the TPC model the fa is typically 
0.9=L0.1, with the small a solutions trending to fa > 2. In contrast, the few acceptable OG models 
have fa < 0.1. The resulting distance for the TPC model is typically d ~ 250 pc (but in some 
cases can be d < 170 pc). For the OG model we expect d m 750 pc. Both cases are small enough 
to be compatible with the small X-ray absorption discussed in Section 13.21 In summary, the TPC 
model has more acceptable solutions, but would imply a very small distance. The relatively small 
parameter space of acceptable OG solutions, on the other hand, is offset by a larger source distance 
and hence a larger Galaxy volume in which such a pulsar could be found. 

In general, the results o f our X-ray analyses are in broad agreement with previous investi- 



gations (jHalpern et al J 120 02), which were based on a factor > 3 smaller photon statistics. Our 
analysis clearly shows that the spectrum of the candidate counterpart is indeed consistent with the 
one of a nearby, thermally-emitting, middle-aged isolated neutron star. 

The best prospect for refining our understanding of the emission from PSR J1836+5925 would 
clearly come from an accurate distance measurement. This seems difficult to obtain, although im- 
proved X-ray spectral measurements and models could help. Alternatively, if 7-ray pulsar spectral 
models can be developed sufficiently, we may be able to connect the softer off-peak spectrum with 
a particular magnetospheric location. In either case, the low power, large characteristic age and 
relatively close dist ance imply that J1 836+5925 is the harbinger of a large population of old and 



weak 7-ray pulsars (jAbdo et al.ll2009al ld) 
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Table 1: Measured and derived timing parameters for PSR J1836+5925 



Parameter Value a 

MJD range 54647.4-55013.0 

Epoch (MJD) 54800 

R.A. (J2000) 18:36:13.75(3) 

Dec. (J2000) +59:25:30.3(6) 

v (Hz) 5.7715514983(4) 

i) (Hz s" 1 ) -4.97(2) xl0~ 14 

Rms timing residual (ms) 1.3 

Characteristic age, r c (kyr) 1840 

^(ergs^ 1 ) l.lxlO 34 
Surface magnetic dipole field strength (gauss) 5-lxlO 11 

7-ray peak separation (A) 0.51 ± 0.01 



"The numbers in parentheses are the la uncertainties derived from the timing model (see Section 2.1). 
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Table 2: LAT 7-ray spectral results PSR J1836+5925 



Phase region 



Parameter 



MJD range 



Value 



54682.7-55013.0 



Phase-averaged 
< d>< 1 



First Peak (FP) 
0.105 < d>< 0.405 



Second Peak (SP) 
0.632 < 6 < 0.904 



Bridge (BR) 
0.459 < 6 < 0.597 



Off-Peak (OP) 
0.053 >(j>> 0.938 



Photon Flux, F100 (ph cm" 2 s" 1 ) 
Energy Flux, G100 (erg cm" 2 s" 1 ) 
Photon Index, V 
Cutoff Energy, E cut0 H (GeV) 
Normalization, K (ph cm" 2 s" 1 MeV" 1 ) 

Photon Flux, F100 (ph cm" 2 s" 1 ) 
Energy Flux, G100 (erg cm" 2 s" 1 ) 
Photon Index, V 
Cutoff Energy, E cutoB (GeV) 
Normalization, K (ph cm" 2 s" 1 MeV" 1 ) 



Photon Flux, F100 (ph cm" 2 s" 1 ) 
Energy Flux, G100 (erg cm" 2 s" 1 ) 
Photon Index, Y 
Cutoff Energy, E cutoff (GeV) 
Normalization, K (ph cm" 2 s" 1 MeV" 



- 1 ) 



Photon Flux, F100 (ph cm" 2 e 
Energy Flux, G100 (erg cm" 2 
Photon Index, V 
Cutoff Energy, £ cu toff (GeV) 
Normalization, K (ph cm" 2 s" 1 MeV" 1 ) 

Photon Flux, F100 (ph cm" 2 s" 1 ) 
Energy Flux, G100 (erg cm" 2 s" 1 ) 
Photon Index, V 
Cutoff Energy, E cutoS (GeV) 
Normalization, K (ph cm" 2 s" 1 MeV" 1 ) 



6.24 ±0.12 
5.91 ±0.08 
.31 ±0.03 
1.27 ±0.11 
1.82 ±0.06 

8.32 ±0.24 
7.97 ±0.17 
.31 ±0.05 
!.31 ±0.17 
2.42 ±0.12 

6.44 ±0.22 
6.40 ±0.16 
.24 ± 0.05 
!.18±0.18 
2.03 ±0.12 

4.47 ±0.28 
4.09 ±0.17 
..17±0.11 
.64 ±0.23 
1.66 ±0.20 



x 10 
x 10 



x 



-7 
10 



x 10" 10 

x 10" 7 
x 10" 10 



x 10 



-10 



x 10" 7 
x 10" 10 



x 10 



-10 



x 10" 7 
x 10" 10 



10" 10 



4.87 ±0.38) x 10- 
3.63 ±0.19) x 10 

1.59 ± 0.11 

2.65 ±0.57 

(9.90 ± 1.2) x 10" 11 



"See the top panel in Figure 1 for a visual representation of the various phase regions. 

'All errors quoted are statistical. In addition, systematic errors of <5Fioo = (± 30%, -10%). #Giq q = (±20%, -10%) 
5Y = (±0.3, -0.1), and SE cutoB = (±20%, -10%) must be taken into account jAbdo et alj|2009lh . 
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Table 3: Radio and X-ray results for PSR J1836+5925 



Wavelength 


Parameter 


Value 


Radio 


Radio flux density at 350 MHz, So.3s{n2y) 


< 55 




Radio flux density at 0.8 GHz, So.sO^Jy) 


< 7 




Radio flux density at 1.4 GHz, Si. 4 (/uJy) a 


< 3 


X-ray 


(0.2-5 keV) Total observed b X-ray flux (erg cm~ 2 s _1 ) 


5.5xl0~ 14 




(0.2-5 keV) Unabsorbed non-thermal flux (erg cm~ 2 s _1 ) 


3.0xl0~ 14 




X-ray blackbody temperature, kT (eV) 


59±I 7 




X-ray blackbody radius (km) 


(i-5tS:2)do.8 c 




X-ray absorbing column, Nh (cm -2 ) 


< 2.7 x 10 20 




X-ray power law photon index, T 


1.7±0.3 




X-ray power law normalization at 1 keV, Npl (ph cm" s~ 1 


ccV -1 ) (5.7±1.0) x 10" 6 



"There is no measured upper limit at 1.4 GHz. This limit is derived from the observation made at 0.8 GHz (See 
Section 3.1) 

^Total observed and unabsorbed flux coincide since the best fitting column density is Nh=Q- 
c d .s is the distance to PSR J1836+5925 in units of 0.8 kpc. 
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0.8 1.0 1.2 
Pulse Phase 



Fig. 1. — Folded light curves of PSR J1836+5925 with a resolution of 32 phase bins per period. 
Two rotations are shown. The top panel shows all events >100 MeV, along with the different phase 
regions labelled: Off-peak (OP), First Peak (FP), Bridge (BR), and Second Peak (SP) regions. The 
horizontal dashed line represents an estimate of the background due to diffuse emission, illustrating 
the high level of off-peak emission being emitted by the source. The lower four panels show the 
light curves in different energy bands. The darker histogram on the second panel from the top 
shows events with E > 5 GeV. 
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Energy {GeV) Energy (GeV) 



(a) Off-peak (OP) (b) First Peak (FP) 




Energy (GeV) Energy (GeV) 



(c) Bridge (BR) (d) Second Peak (SP) 

Fig. 2. — Energy spectra of the four identified phase regions of PSR J1836+5925. The data points 
represent the measured fluxes obtained from likelihood fits in different representative energy bands 
where the pulsar is modeled as a power law, while the line shows the best-fit model obtained in the 
unbinned maximum likelihood analysis over the entire energy range, along with the la "bowtie" 
confidence region. Top Left — Off-peak (OP). Top Right — First Peak (FP). Bottom Left — 
Bridge (BR). Bottom Right - Second Peak (SP). 
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Phase 

Fig. 3. — Cutoff energy (top) and photon index (bottom) as a function of phase for PSR J1836+5925 
using 15 equal-count bins containing ~650 events each. The dashed line in both panels shows the 
>100 MeV folded light curve of PSR J1836+5925. 
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Fig. 4. — Flux of PSR J 1836+5925 as a function of time in 5-day time bins, showing no evidence 
for variability (see Section 2.4.2). 



-24- 



O r 



> CO 

O 

i 

m 

E o 

o i- 

> 

CD 

°? 
O 



i 

O 



i — i — i i 1 1 1 ii 1 — i — i i ■ ■ i ■ i 1 — i — i i 1 1 1 ii 1 — i — i i 1 1 1 ii 1 — i — i i 1 1 mi 1 — i — i i 1 1 mi 1 — i — i ■ ■ i ■ 1 1 1 — i — i i 1 1 1 ii 1 — i — i i 1 1 1 ii 




PSRJ1 836+5925 



_l I ' I I I I I I I I ' I I I I 



1 000 1 4 1 o 5 
Energy (keV) 

Fig. 5. — Unfolded XMM-Newton spectrum of PSR J1836+5925, compared to that of Geminga, 
whose data have been reprocessed to take advantage of the new calibration files which cover the 
energy range down to 0.15 keV (EPIC status of calibration and data analysis Document XMM-SOC- 
CAL-TN-0018 |http://xmm2.esac.esa.int/docs/documents/CAL-TN-001 8.pdf). The best fit models 
are superimposed (for a thorough report on Geminga, see lCaraveo et al.ll2003l ). We extrapolate the 
models out to the LAT energy range and show the best-fit phase-averaged LAT spectra for the two 
pulsars. The overall similarity is apparent. 



